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Abstract. 


Stable baryonic Q-balls, which appear in supersymmetric extensions of the 
Standard Model, could form at the end of cosmological inflation from fragmentation of 
the Affleck — Dine condensate. We reconsider astrophysical constraints on such Q-balls 
as dark matter candidates. Baryonic Q-balls interact with matter by absorbing the 
baryon number and, effectively, leading to a rapid baryon number non-conservation. 
We have recently shown that this process can occur at a much faster rate than that 
used in previous calculations. As a consequence, stability of neutron stars imposes a 
stringent constraint on the types of Q-balls that can be dark matter. Only the Q-balls 
that correspond to baryonic flat directions lifted by baryon-number violating operators 
are allowed as dark-matter candidates. 
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1. Introduction 


Supersymmetry (SUSY) provides two types of dark matter candidates: the lightest 
supersymmetric particle (e.g., neutralino or gravitino), and stable Q-balls. SUSY Q- 
balls are nontopological solitons that carry baryon number B]. They are either 
stable or have lifetimes in excess of the age of the universe in theories with gauge- 
mediated SUSY breaking [4]. At the end of inflation, large SUSY Q-balls can form via 
fragmentation of the Affleck-Dine condensate [5] [6] [7]. 

SUSY Q-balls can interact with matter fermions. A quark scattering off a baryonic 
Q-ball can convert into an antiquark, while increasing the baryon number of the Q-ball. 
Hence, Q-balls catalyze baryon number non-conservation in the sector of ordinary, non- 
supersymmetric particles. The effective B-nonconservation catalyzed by Q-balls was 
considered in Ref. [S] in connection with astrophysical limits on Q-ball dark matter. 
However, it was recently shown [9] that the rate of this process is considerably higher 
than that used in Ref. [8]. In this paper we reconsider astrophysical limits on dark 
matter Q-balls and show that a wide variety of them are ruled out because they would 
destabilize neutron stars on cosmologically short time scales. 


2. Interactions of SUSY Q-balls with matter 


In a previous paper we reanalyzed the interaction of Q-balls with ordinary matter. 
We found that quarks falling on a Q-ball are reflected as antiquarks with a probability 
of order one, practically independent of the parameters of the theory. In other words, 
Q-balls convert the matter into antimatter on their surface (or antimatter into matter, if 
placed in an anti-matter environment). Baryon number is conserved during this process: 
after reflection of an antinucleon, the baryonic charge of the Q-ball increases by 2 units. 

Having described the interaction of SUSY Q-balls with nucleons, we now reexamine 
some bounds on dark-matter Q-balls. The laboratory limits do not change significantly 
from the earlier estimates [8]. This is because nucleons reflected from the Q-ball surface 
as antinucleons have a very small mean free path in the surrounding matter. They 
annihilate quickly, mainly into pions. The energy release is somewhat higher than it 
was thought to be [S|, but this is an order one effect, which does not make a difference 
for most bounds. 

However, the analysis of astrophysical bounds, in particular those from Q-balls 
captured in the neutron stars, must be done differently. The rate of baryon number 
absorption by a Q-ball was underestimated in Ref. [I0] by several orders of magnitude. 
One can expect a dramatic change in the astrophysical bounds because of this. In fact, 
such rapid absorption of nucleons requires that we reexamine the inner fluid dynamics 
of the neutron star. Although the dynamics of nuclear matter is difficult to describe, 
we understand it well enough to know that certain types of Q-balls would consume a 
neutron star in a few thousand years. This is far less than the lower bound lifetime 
of neutron stars and suggests that such Q-balls must not exist in large amounts in the 
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universe. 

We will also address another important issue overlooked in the earlier work. Non- 
renormalizable operators alter the flat directions in the potential and can modify the 
shapes of solutions for SUSY Q-balls. We will consider interactions of Q-balls of different 


types. 


3. SUSY Q-balls and the flat directions 


Stable Q-balls develop a VEV along some flat direction. This flat direction can be 
parameterized by a scalar field y. The VEVs of squarks and sleptons are proportional 
to y. In theories with gauge-mediated breaking of supersymmetry, the effective potential 
along this flat direction is nearly constant, V (p) x Af, up to some Ymax Which depends 
on the type of the flat direction +. There exists a classical spherically symmetric solution 
of the field equations, a non-topological soliton, in which 


p= dof (r), (1) 
where f(r) ~ sin(wr)/(wr) for r < R = m/w and zero for r > R. It carries the global 
charge Q. If ġo < Ymax, different parameters depend on Q as follows: 


po ~ AQ, M~ AQF, wr AQ’. (2) 


We shall call the non-topological solitons of this type flat direction (FD) Q-balls. 

Flat directions are in general lifted by higher dimension non-renormalizable 
operators and the effective potential starts to grow at y > Ymax. Thus, when 
the charge Q exceeds some critical value Qe, the properties of Q-balls change. The 
details of this analysis will appear elsewhere [12], here we reproduce the essentials only. 

Generally, the lifting terms can be written in the form 


N é od n—-l+m o“ n—l-m 
V (P)iting S An M (£) (5) , (3) 


where the smallest possible n and corresponding to it m are some integers that are 
determined by the structure of the flat direction and can be found from Tables 4 and 
5 of ref. [O]. We will assume that An is of the order one and that M is the grand 
unification scale ~ 10'6GeV and that A ~ 10° GeV. 

The relations (2) are valid only if the lifting potential is small compared with A‘*. 
So, the Q-ball properties must be re-analyzed when their charge reaches the critical 


value Q. determined from the condition V"(|])iiting ~ A4, and equal to 


4n—12 


Qaa (DT), (1) 


It can be seen from [II] that most flat directions are lifted by monomials of dimension 


n = 4. At this low suppression, only Q-balls with small charge Q < Qe ~ 10!" can 
retain their « Q°/4 energy dependence despite the lifting. On the other end, all flat 


{+ We omit possible log factors, because they are not essential for the discussion. 
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directions are lifted by monomials of dimension n = 7, which corresponds to a critical 
charge Q. ~ 10%. 

The fate of Q-balls with charges larger than Qe depends crucially on the value of 
the integer m in (B). For m Æ 0 the lifting potential not only breaks the degeneracy 
of the flat direction, but introduces an explicit breaking of the global U(1) symmetry 
responsible for Q-ball existence. For these flat directions non-renormalizable operators 
induce strong baryon number non-conservation, which leads to fast Q-ball decay. As 
the value of Qe happens to be much smaller than the baryon number of a star, no 
astrophysical limits on this type of Q-balls can be derived. These Q-balls absorb only a 
small fraction of the baryon material of a star and disappear afterwords. In notations 
of ref. [I] Q-balls related to flat directions QLe, QLd, Lude, QLde, QLud, QLude 
are of this type. 

If, on the contrary, m = 0, the lifting term does not break the global U(1) symmetry 
and the Q-ball charge can grow beyond Qe, but the dependence of the parameters of 
the Q-ball on Q is modified. Now, with the lifting terms incorporated, the quantity 
V(y)/|y|? does have a minimum at non-zero finite value of y, and the analysis of [I] is 
applicable. Thus, 


Po ~ Pmax M ~ uQ, wrp, Rr (QAP, (5) 
where u’ ~ V(Qmax)/Priax: Changing n from n = 4 to n = 7 the frequency w ranges 
from 10 MeV to 107? MeV and the VEV ġo from 10’ to 10! GeV. In notations of 
ref. [T] Q-balls related to ud, ue, QL, ude flat directions are of this type. In the 
following we will be interested in these Q-balls only. 

Though in general the critical charge Qe depends on the structure of the flat 
direction and on parameters of the underlying unified model [12], the exact values of the 
critical charge and of w are of little importance for us. What is essential is that (i) Qe 
is much smaller than the baryon charge of the neutron stars Qas ~ 10°"; (ii) w is much 
smaller than the nucleon mass and smaller than the chemical potential for the baryon 
number in the center of a neutron star; and (iii) that ġo > w for dark matter Q-balls. 
These types of non-topological solitons shall be dubbed CD (curved directions) Q-balls. 

Very large Q-balls can become black holes. For FD Q-balls this happens when the 
global charge reaches 


TeV \* 
Qin = 8-10 (L), (6) 
A 
and for CD Q-balls the corresponding number is 
TeV\? (M 
a 22-108 (EE) (8) ; 
A w 


While for FD Q-balls it exceeds the typical baryonic charge of a star, for CD Q-balls it 
is less than Qys, which is important for discussion below. 
A summary of the Q-ball properties can be found in Table I. 
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Table 1. Summary of properties of Q-balls of different types 


4. Astrophysical limits 


We now ask whether such a large processing rate of baryonic matter by the relic SUSY Q- 
balls could affect some astrophysical observations which could be used to either detect 
dark-matter Q-balls or rule them out. Q-balls are extremely dense and heavy and, 
therefore, they do not stop in most astrophysical objects, such as stars or planets [IQ]. 
They would stop, however, in both neutron stars and white dwarfs. Furthermore, simple 
estimates suggest that most neutron stars and white dwarfs are infected by Q-balls 
within a short period of time [10]. Neutron stars and white dwarfs are, therefore, the 
most reasonable places to look for Q-ball signals. 

A Q-ball captured by a neutron star or a white dwarf begins consuming baryonic 
matter. At some pooint, the neutron star mass can go below its stability limit, which 
is about 0.2 Mo [I3]. Then the neutron star explodes. This happens when a CD Q-ball 
inside becomes a black hole. If the time scale for this process is less than the known 
ages of some neutron stars and, then such Q-balls must not exist. Spindown rates of 
some pulsars is as low as (P/P) ~ (0.3 — 3) x 107!yr7!, which corresponds to the age 
in billions of years. Some pulsars are also known to be (at least) as old as 10 Gyr based 
on the cooling ages of their white dwarf companions [14]. One must, therefore, rule out 
all Q-balls that would destroy neutrons stars in less than 10!° years. 

As we have already mentioned, the estimate of the neutron star lifetime made 
in Ref. [10] was based on an incorrect rate of the baryon number absorption by non- 
topological solitons. In this section we will revise and correct this analysis. We shall 
also incorporate two other effects that were not essential for the analysis in Ref. [I0] 
because of the assumed small rate of the baryon number transfer from the neutron 
matter to Q-balls. The first effect is related to the fact that the rate of baryon number 
transformation on the surface of Q-balls is large, and, therefore, there is a large pressure 
from the radiated pions. In this case the rate of absorption is determined by the fluid 
dynamics and relative pressures of neutrons and pions rather than by the probability of 
matter-antimatter transformation on the Q-ball surface. The second effect is that the 
charge-energy relations are changing after the Q-ball reaches the critical charge discussed 
in Section 2, see eq. (6). Unfortunately, relatively little is known about the internal 
structure of neutron stars and white dwarfs. Therefore, we will present two simplified 
models to get an idea of the timescale for neutron star and white dwarf consumption. 
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4.1. Unitarity bound 


Let us present a lower bound on the lifetime of a neutron star infested with Q-balls. 
Let us assume that the baryon number density near the Q-ball surface is given by the 
nuclear density no ~ 4- 1073 GeV? at the center of the neutron star without a Q-ball 
[3]. Then the absorption rate simply coincides with the number of baryons falling on 
the Q-ball, 


= TTO NU ; (8) 


dt 
where v is the typical nucleon velocity, taken to be of the order of the speed of light in 


what follows. 

Let us now assume that the nuclear density in the center of the star does not change 
with time (this is a conservative assumption because the central density decreases when 
the baryon number of the neutron star is consumed by a Q-ball). Then the lifetime of 
a star populated by FD Q-balls is of the order of 

2 
t~ A (Qua)? ~ 10° years . (9) 
No 
For the CD Q-balls (a more realistic case, as the FD Q-balls are transformed into CD 
Q-balls after their charge increases past Qe) the lifetime is much smaller, 
2 42/2 
pm Cees Oua) ~ O(1) day (10) 

In contrast, a white dwarf is about 10° times less dense than a neutron star, so the 
time scale for a white dwarf to be consumed is about 10" years for a FD Q-ball and 
about 10° years for a CD Q-ball. 

All but one of these estimates are considerably smaller than the typical ages of 
neutron stars or white dwarfs. We will see, however, that the baryon number absorption 
is mainly determined by the transport properties of quarks and nucleons in the interior 
of the neutron stars and that the speed of baryon number consumption by Q-balls is 
in fact considerably smaller. Nevertheless, it happens to be large enough to completely 
exclude Q-balls (with m = 0 in eq. (8)) as candidates for dark matter. 


4.2. Hydrodynamic Considerations 


It is clear that the unitarity limit is unreasonable. However, it is not clear what the 
appropriate transport mechanism inside the neutron star will be. We know that near 
the Q-ball neutrons will be consumed and turned into pions (either directly at the Q- 
ball surface, or in an annihilation). We also know that a neutron star is held stable 
by the interplay of gravity pulling inward and the tremendous degeneracy pressure of 
the neutrons pressing outward. If the change in neutron density results in a pressure 
change, then the transport will be dominated by pressure driven hydrodynamic flows. If 
instead, the pressure remains constant, then the transport will be dominated by either 
convection or diffusion. Therefore, we must first determine whether or not the pressure 
can remain constant near the Q-ball. 
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There are two simple ways to estimate the pressure near the center of the star. 
Macroscopically, we can determine the required pressure to maintain the star in 


hydrostatic equilibrium. For a constant density the pressure is 
3 M? 1 
P=— —. 11 
87 Mp? Rt (E 
Using M ~ 10°” GeV, and R = 10km ~ 5x10!GeV~ yields 2x1074GeV* ~ (.12Gev)*. 
In the second approximation we use the pressure from a degenerate non-relativistic 


fermi gas. 
2)1/3 
(67°) në. 
5m 
Using the density of the star as 1015 g/cm we find that P œ 10-4GeV* ~ (.1GeV)4. 
While both approximations are highly simplified, they give the same order of 


(12) 


magnitude for the pressure which we assume is reliable. As the Q-ball consumes baryon 
number, the degeneracy of neutrons will no longer be able to maintain this high pressure. 
The pressure will then have to be maintained either thermally, or by pions before they 
are able to decay into other particles. We will see shortly that neither option is viable. 


4.2.1. Thermal Pressure. The thermal pressure from light particles is 
4 


Px I (13) 
where g is the number of light degrees of freedom. In our case these are 2 photons, 4 
electrons, 2 neutrinos and possibly 3 pions (depending on the exact temperature). This 
suggests that we will have an extrememly high temperature of about 100 MeV. 
Such high temperatures can not be maintained because neutrinos carry the energy 
away faster than the Q-ball can produce it by absorbing neutrons. Therefore we know 
that the pressure near the Q-ball can not be maintained by thermal equilibrium. 


4.2.2. Pion Pressure. While the pressure can not be maintained by thermal effects, 
we may hope that it is maintained by the extra pions produced before they are able to 
decay. Because the pions interact on strong length scales, but decay over electromagnetic 
and weak time scales, there is a significant time when they behave much like thermal 
conserved particles. 

We will therefore assume that near the Q-ball, we have a mixture of neutrons and 
pions, each with an average kinetic energy of about .1 GeV, and that the pressure of 
either component is simply the number density times the kinetic energy. The pions 
will continuously change between species, but since the average kinetic energy is lower 
than their mass, no new pions will be produced. We shall also assume that once a 
pion decays, it no longer plays a significant role in the pressure of the system. This 
is because photons, neutrinos, and electrons all have much smaller cross sections than 
neutrons and pions, and they will immediately leave the area near the Q-ball. Any 
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such particles entering the Q-ball region will have the much lower energies typical of the 
colder outside region. 

In this situation, the pions will cluster near the Q-ball where they are produced, 
but they will wander away from the Q-ball essentially due to random walk diffusion. 
However, the distance they reach is limited by their lifetime 7. Thus, they will reach an 
average distance from the Q-ball of about /3AT We will model their density profile as a 
falling Gaussian with this decay length. A one dimensional model is sufficient because 
the decay length is much less than the radius of the Q-ball. 


Nr (£) = Npi(0) exp (-=) l (14) 


Then the total rate of pion loss is, approximately, 


dN, = 1 2 [2 r2 
T = rro f nO) exp ( 7 Ja 


T 


= 27/4 | A ng(0)rg? (15) 


To maintain pressure, n,(0) needs to be about the same as the density of neutrons 


far from the Q-ball, but still near the center of the star. That is no ~ .1GeV’. We 
-1/3 
0 . 


will also take \ ® n That is we are letting the mean free path be about the 
spacing between particles. Finally r œ~ 10°GeV~'. (This is the lifetime of a neutral 
pion. We assume that the decay rate is dominated by the electro-magnetic decays of 
neutral pions.) From this we find 
dN, 
dt 


This is the rate at which pions are lost to decays. (If we require them to be lost 


= 107% Q?/3GeV . (16) 


through weak decays, the number changes to about 107%. This is less of a change than 


-1/2 not as 7: 


we might expect because the lifetime only enters as T In essence, the 
pions are decaying slower, but there are more of them. However, for such a situation to 
happen we must justify a lack of neutral pions. Given the frequent strong interactions 
this seems unlikely.) Now, one neutron can only supply about 4-5 pions, so the rate of 
baryon consumption required to maintain pressure near the Q-ball is no more than one 
order of magnitude below this. That means 

dNn 


dt 
t = 3 x 10" Qpgy! GeV} ~ 1500 years . (18) 


~ 1071 Q? GeV , (17) 


This is clearly much longer than the unitarity limit, but still well below the known 
lifetimes of neutron stars. Because they are smaller, grow more slowly, and must 
consume a larger number of neutrons, the lifetime of a FD Q-ball is much longer, about 
10!° years. 

The neutron star lifetimes are summarized in Table IV. 
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| | FD Q-balls | CD Q-balls 
1500 years 


Table 2. Neutron star lifetimes for Q-balls of different types. 


Of course in reality, the neutron density near the Q-ball will be decreased slightly 
due to their rapid consumption. However, as we have shown this will result in a drop 
in pressure which will cause more neutrons to flow to the site. If the density dropped 
by several orders of magnitude, this would greatly increase the lifetime of the star. 
However, the pressure is so great that a slight (order 1 or less) drop in neutron density 
near the Q-ball creates a pressure gradient that is sufficiently large to replenish the 
consumed neutrons. Such a drop would not greatly change the lifetime of the neutron 
star. 

The predicted lifetime of CD Q-balls is so small that we can clearly rule out their 
existence. If they existed there would be no neutron stars older than a few thousand 
years. FD Q-balls could still survive the astrophysics bound. However, as mentioned 
earlier, large FD Q-balls are unlikely to exist since all the flat directions are lifted for 
VEVs much smaller than those for which the Q-ball reaches the size of a neutron star. 

However, this limit applies only to Q-balls with m = 0, which conserve baryon 
number at all values of Q. Q-balls for which m 4 0 do not grow beyond a certain size 
because the baryon number non-conserving interactions in their interior prevent Q from 
growing beyond the critical value Qe. Such Q-balls remain small, and, therefore, they 
do not consume a neutron star on cosmologically interesting time scales. 

We note in passing that the corresponding analysis for the white dwarfs would 
be more complicated because the main pressure supporting a white dwarf against 
gravity comes from the Fermi pressure of degenerate electrons. This pressure is coupled 
only electromagnetically to the strong dynamics of nuclei and Q-balls. The surface 
temperature of the white dwarfs could be used to set a limit on the amount of heat 
generated inside by the Q-ball [I5]. There are many 10 Gyr old white dwarfs that have 
cooled to very low temperatures and have luminosities as low as 


Lwa = 3 x 10% Lo = 7 x 10” erg/s. (19) 


However, we do not need to consider white dwarfs in detail because, for m = 0, Q- 
balls are already ruled out by the stability of the neutron stars. For m # 0, the rate 
of energy release from consumption of nuclear matter is much lower than the value in 


equation (I9). 


5. Conclusions 


We have revised the bounds on the relic dark-matter Q-balls based on the existence 
of neutron stars and their stability on the time scale of billions of years. Assuming 
SUSY Q-balls have a sufficient abundance to be the dark matter, relic SUSY Q-balls 
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are captured by a neutron star shortly after its formation. The Q-ball consumes ambient 
neutrons and adds to the baryon number of its condensate, while dissipating the excess 
energy in neutrinos and photons. Meanwhile, the neutron star mass decreases. For most 
SUSY Q-balls, this process is too rapid, and it can destabilize a neutron star on time 
scales much shorter than a billion years. Hence, such Q-balls are ruled out. The only 
SUSY Q-balls that survive this bound and remain dark-matter candidates are those 
that have their flat directions lifted by the baryon number violating operators. Such 
Q-balls cannot grow beyond a certain size, which turns out to be too small for a rapid 
consumption of nuclear matter. 
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